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75.85. +t - Multiferroics 

75.60. Jk - Magnetization reversal 

75.78. -n - Magnetization dynamics 


t - Electric field induced magnetization switching in single-phase multiferroic materials is 
intriguing for both fundamental studies and potential technological applications. Here we develop 
a framework to study the switching dynamics of coupled polarization and magnetization in such 
multiferroic materials. With the coupling term between the polarization and magnetization as 
an invariant dictated by the Dzyaloshinsky-Moriya vector, the dynamical systems study reveals 
switching failures and oscillatory mode of magnetization if the polarization and magnetization 
relax slowly during switching. 
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Introduction. — Multiferroics usually represent ma¬ 
terials that are both ferroelectric and ferromagnetic (THE]. 
Such materials in single-phase were usually thought to be 
rare [13], and hence multiferroic composites in 2-phase, 
i.e., a ferroelectric layer strain-coupled to a ferromagnet, 
are usually deemed to be the replacement |5lfT7l20] . How¬ 
ever, there have been recent resurgence of interests EZXE2I 
and some mechanisms of coupling polarization and magne¬ 
tization in single-phase materials are coming along [23)425] . 
This can lead to possible technological applications [26] of 
switching a bit of information (stored in the magnetiza¬ 
tion direction) by an electric field [77] . This eliminates 
the need to switch magnetization by a cumbersome mag¬ 
netic field or spin-polarized current [13], although new con¬ 
cepts are being investigated e.g., utilizing giant spin-Hall 
effect [28] . The electric field switches the polarization and 
the intrinsic coupling between the polarization and mag¬ 
netization switches the magnetization between its 180° 
symmetry equivalent states. One way to couple polar¬ 
ization and magnetization that has taken attention is due 
to Dzyaloshinsky-Moriya (DM) interaction [29j[30], which 
arises due to spin-orbit correction to Anderson’s superex¬ 
change m- This is called ferroelectrically induced weak 
Ferromagnetism (wFM), in which two magnetic sublat¬ 
tices of an antiferromagnet cant in a way to produce a 
residual magnetization [ 231137435] . 
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While first-principles calculations and experiments have 
been underway on the search of strongly-coupled multifer¬ 
roic magnetoelectric materials possibly working at room- 
temperature, little have been studied on the dynamical 
nature of switching. The study of switching dynamics of 
magnetization in multiferroic composites, i.e., a piezoelec¬ 
tric layer strain-coupled to a magnetostrictive nanomag¬ 
net, have been very successful to understand the perfor¬ 
mance metrics, e.g., switching delay, energy dissipation, 
and switching failures mmm Here, the switching dy¬ 
namics of polarization is studied by forming a Hamilto¬ 
nian system with Landau-Ginzburg functional [36], while 
the magnetization dynamics is studied by the Landau- 
Lifshitz-Gilbert (LLG) equation of motion [37]EE]- We 
focus on magnetization switching due to electric field in¬ 
duced polarization switching, i.e., converse magnetoelec¬ 
tric (ME) effect for technological applications rather than 
the switching dynamics due to direct ME effect. We par¬ 
ticularly consider the dynamics in single-domains with an 
eye to achieve high-density of devices rather than consider¬ 
ing domain walls in higher dimensions, for which we need 
to consider the competition between the exchange interac¬ 
tion and dipole coupling among the spins [35] . Note that 
switching dynamics in BiFeOs has been studied using a 
first-principles-based effective Hamiltonian within molec¬ 
ular dynamics simulations mm- Here we perform a com¬ 
prehensive analysis in emerging strongly-coupled multifer¬ 
roics (ferroelectrically-induced wFM by design dictated by 
DM interaction [23]) by varying different parameters in 
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Fig. 1: (a) Energy of the two spins in a predominantly antiferromagnetic configuration with respect to canting angle 0 C . The 
canting happens due to Dzyaloshinsky-Moriya (DM) interaction. The DM vector is proportional to polarization, i.e., D oc P 
and acts in the ^-direction. When 0 C is positive, the net magnetization M points up (+y- axis), and if 0 C is negative, the 
net magnetization M points down (—y-axis). Without any canting, the net magnetization is zero as in the case of a perfect 
antiferromagnet. (b) Magnetization initially points along the —y-xis. To respect the invariant dictated by the DM interaction 
[oc P • (L x M)], with the reversal of polarization P, either the magnetization M (case 1) or the AFM vector L (case 2) may 
flip. 


the LLG dynamics. The analysis of switching dynamics 
reveals very significant motion of magnetization when po¬ 
larization is switched by an electric field. It is shown that 
magnetization may fail to switch or even can go into an 
oscillatory state of motion. The phenomenological damp¬ 
ing parameter for both polarization and magnetization 
plays a crucial role in shaping the dynamics of the coupled 
polarization-magnetization in these multiferroic materials. 

Model. — We consider two spins, one representative 
to each magnetic sublattice of an antiferromagnet, to build 
up the present model. The dynamics of the two spins Si 
and S 2 can be described by the Landau-Lifshitz-Gilbert 
(LLG) equation [371138] as follows: 


K is the single-ion anisotropy constant, 7 ' = 7/(1 + <a 2 ), 7 
is the gyromagnetic ratio of electrons, a is the phenomeno¬ 
logical Gilbert damping constant [38], and S = |Si| = 

| S 2 1. As required, it is possible to include the long-range 
interaction too in the energy term m- The net magneti¬ 
zation M and the antiferromagnetic (AFM) vector L for 
the two spin system are M = Si + S 2 and L = Si — S 2 , 
respectively. Note that the following two identities hold: 
M • L = 0 and M 2 + L 2 = 4S 2 . 

The polarization dynamics is based on the Landau- 
Ginzburg functional [36] 


G = 


p2 + ^2 p 4 

2 4 


E.P, 


(5) 


^ = -| 7 , |S 1 xHsi-^S 1 x(S 1 xH S i) (1) 

^ = -|7 , |S 2 xH S 2-^S 2 x(S 2 xHs 2 ), ( 2 ) 

where Hsi and h S2 are the effective fields on the spins Si 
and S 2 , respectively, defined as Hsi = —{dH/d Si) and 
Hs2 = —{d'Hjd S 2 ), T~L is the potential energy of the two 
spin system, expressed as 

n = -JSi • S 2 — D • (Si x S 2 ) - KS'l z - KSl z , (3) 

J denotes the exchange coupling between the spins, D is 
the Dzyaloshinsky-Moriya (DM) vector [291130] (here, we 
will consider the case when the vector D points along per¬ 
pendicular to the plane (x-y plane) on which the spins 
reside, i.e., along the ^-direction and proportional to po¬ 
larization P, which is also in the ^-direction [23]) expressed 
as D = D(t) e z making 

D(t) oc P(t), (4) 


where a\ and a 2 are the ferroelectric coefficients (both are 
greater than zero) and E = E e z is the applied electric field 
that switches the polarization P in the ^-direction. We as¬ 
sume single-domain case [13] and follow the prescription 
in Ref. [36] to trace the trajectory of polarization. Note 
that polarization is switched by moving ions, which cou¬ 
ples to the magnetization dynamics via the DM term D 
[see Eq. ©]• On the other hand, rotation of spins does 
not quite move the heavy ions affecting the polarization. 

Figure [lji depicts how the asymmetric Dzyaloshinsky- 
Moriya (DM) interaction can lead to two anti-parallel 
magnetization directions, i.e., 180° symmetry equivalent 
states. Depending on the sign of the canting angle 0 C of 
the spins, the direction of the DM vector changes and 
the energy expression as in the Equation © gives rise to 
two magnetization states in opposite directions. The DM 
vector D is proportional to polarization P and hence, if 
we switch the polarization, two cases can happen to re¬ 
spect the invariant due to DM interaction P • (L x M) 
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Polarization (nC/cm 2 ) time (ns) 

Fig. 2: (a) The potential landscape of polarization with electric field as a parameter [Equation JJJ]. Note that it requires a critical 
electric field to topple the barrier between polarization’s two 180° symmetry equivalent states, (b) Switching of polarization 
with the application of an electric field. Initially, the polarization was pointing towards —z- axis. With the application of electric 
field, the polarization does not reach instantly towards the +z-axis, how fast the polarization relaxes to the minimum energy 
position depends on the polarization damping. Note that after the withdrawal of electric field, the polarization direction is 
maintained, i.e., the switching is non-volatile. The slight increase of polarization over P s is due to the application of electric 
field [see the potential landscapes in part (a)], which can be followed from the Equation (|5} too. 


[oc D • (Si x S 2 ) term in Equation (|3J)]: (1) The magneti¬ 
zation M can change the direction (i.e., switches success¬ 
fully), and (2) The AFM vector L may change the direc¬ 
tion (i.e., M fails to switch). The two cases are depicted 
in the Fig. [lj). 

Results and Discussions. — We consider a per- 
ovskite system NiTiOs mm in R3c space group 05] as 
a prototype to analyze the switching dynamics. Although, 
NiTiOs in R3c space group is not yet experimentally re¬ 
alized, the concept of polarization-magnetization coupling 
predicted from group theory is promising. The parame¬ 
ters are chosen as follows: saturation polarization P s = 
110 /jl C/cm 2 , ferroelectric coefficients a\ = 1.568 x 10 10 
Vm/C, Ci 2 = 1.296 x 10 10 Vm/C, polarization damp¬ 
ing p = 0.286 VmSec/C (that switches the polarization 
in realistic time 100 ps [46], see Fig. E h S m 1.6/is, 
M = 0.25 /is, J = —2.2 meV, D s = 0.35 meV [corre¬ 
sponding to P s , i.e., D(t) = (D s /P s ) P(t)\, K = —0.03 
meV [23] . We will consider that the electric field switches 
the polarization from — P s to -\-P s in the ^-direction. 

The magnetization damping, through which magne¬ 
tization relaxes to the minimum energy position, can 
have a wide range of values (10 -4 - 0.8) by modifying 
the spin-orbit strength, doping etc. and it can be de¬ 
termined by ferromagnetic resonance (FMR), magneto¬ 
optical Kerr effect, x-ray absorption spectroscopy, and 
spin-current driven rotation with the addition of a spin- 
torque term 0TH19] . Hence, we focus on investigating the 
magnetization dynamics for a wide range of phenomeno¬ 
logical damping parameter. 

We will initially assume the single-ion anisotropy K = 
0 and we will see later the consequence of considering 
it. Figure [3] shows the dynamics of magnetization when 
damping parameter is on the higher side, e.g., 0.1. We 


see that magnetization has switched successfully in the 
end (see Fig. [3]i), while the AFM vector did not switch 
(see Fig. [3t). Note that the spins Si and S 2 are deflected 
from the x-y plane, in the ^-direction due to rotational 
motion of magnetization. Also, note that magnetization’s 
x- and ^-component and AFM vector’s ^/-component have 
not changed at all due to the complimentary dynamics of 
the spins Si and S 2 . This corresponds to the case (1) in 
Fig. [T 

Figure 0] plots the dynamics when magnetization damp¬ 
ing a = 0.01. We see that magnetization has failed to 
switch (see Fig.®!), while the AFM vector is switched suc¬ 
cessfully (see Fig. m- Magnetization was on the way to 
change its direction, but eventually, magnetization came 
back to its initial state. This corresponds to the case (2) 
in Fig. [TJd. Due to low damping, we notice ringing in all 
the plots in the Fig. 0] 

For the lower damping of a = 0.01, from the simulation 
results as shown in the Fig. 0] the positions of the two 
spins have just got interchanged, which is depicted as the 
case ( 2 ) in Fig. [TJd. Since the canting angle of the spins are 
small (0 C ~ 5°), one can say that the spins have rotated 
much more than that of the case for the higher damping 
of a = 0.1 [case (1) in Fig. [TJd and the simulation results 
as shown in the Fig. [3]. While both the cases as shown in 
the Fig. [TJd respect the DM invariant at steady-state , the 
dynamics of magnetization dictates the final state that is 
reached. With a lower damping, the spins get deflected 
out-of-plane more (see the ^-components of the spins Si 
and S 2 in the Figs. [3] and 0]) and this out-of-plane excursion 
eventually can lead the spins to interchange their positions 
as can be noticed from the Fig. 0] The interchange of the 
spins Si and S 2 indeed respects the DM invariant [case 
(2) of Fig. [1 ]d] , but the magnetization M fails to switch in 
this case, while the AFM vector L gets switched. 
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Fig. 3: Dynamics of magnetization for damping parameter a = 
0.1. Magnetization does switch successfully, (a) Dynamics of 
Si, (b) Dynamics of S2, (c) Dynamics of AFM vector L, and 
(d) Dynamics of magnetization M. 

To understand the magnetization dynamics further that 
how switching may be successful even at magnetization 
damping a = 0.01, we first investigate its dependence on 
the polarization dynamics. Simulation results show that 
magnetization switches successfully if we make the polar¬ 
ization damping 200 times faster [see Fig. [5^]. Basically, 
due to the coupling between the polarization and mag¬ 
netization, if polarization is switched faster, magnetiza¬ 
tion is also switched faster, which makes the switching 
successful. We further investigate the effect of single-ion 
anisotropy parameter K on magnetization dynamics. The 
ion-anisotropy basically adds an extra field that tries to 
keep the magnetization in-plane (i.e., x-y plane) and the 
simulation results show that magnetization switches suc¬ 
cessfully if we take the single-ion-anisotropy into account 
and the ringing in the magnetization dynamics does not 



time (ns) 

Fig. 4: Dynamics of magnetization for damping parameter a = 
0.01. Magnetization fails to switch successfully, (a) Dynamics 
of Si, (b) Dynamics of S2, (c) Dynamics of AFM vector L, 
and (d) Dynamics of magnetization M. 

show up in this case [see Fig. Eb]. However, if the single¬ 
ion-anisotropy is reduced to a value of K — —0.003 /ieV, 
it is noticed that the magnetization fails to switch success¬ 
fully. 

We further study the effect of varying the DM interac¬ 
tion strength D s on the switching dynamics. It is found 
quite obviously that as D s decreases for a fixed canting an¬ 
gle 0 C , i.e., polarization-magnetization coupling weakens, 
the AFM vector L deflects more, which can be interpreted 
as that the magnetization M is more prone to switching 
failures. For a = 0.1, if D s is reduced 10 times, the mag¬ 
netization M still switches successfully. 

An interesting investigation would be to see whether 
magnetization, being a rotational body, oscillates for a 
certain range of damping parameter. For example, a spin- 
polarized current can spawn oscillatory states in a nano- 
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(a) Dynamics of Magnetization M 




time (ns) 


Fig. 5: Magnetization switches successfully at damping parameter a = 0.01. (a) Polarization is switched 200 times faster by 
changing the polarization damping. The magnetization switches successfully, (b) The single-ion anisotropy constant K is taken 
into account. In this case, magnetization switches successfully too. 


magnet [50]. Figure [6^ shows the magnetization dynam¬ 
ics when damping parameter a = 0.001. The single-ion- 
anisotropy is included here and magnetization was able 
to get past towards the +//-direction, however, could not 
settle there and oscillates with a time period of 3.3 ps. 
With further lowering of the damping parameter, magne¬ 
tization still oscillates, however with a lower frequency [see 
Fig. [6)3]. 

The oscillatory mode of magnetization too can be ex¬ 
plained from the out-of-plane excursion of the spins due to 
low damping. The spins get deflected out-of-plane (i.e., z- 
direction) and when they go completely out-of-plane, they 
continue rotating and reach the out-of-plane in the oppo¬ 
site directions than the previous ones. Therefore the spins 
sustain a self-oscillation. The DM interaction ensures that 
the x- and ^-component of the spins are canceled out and 
the //-components are added due to symmetry, as can be 
noticed in the Fig. [6] Note that such self-oscillation oc¬ 
curs and sustains even in the absence of an external electric 
field, making the system unstable. Such spontaneous self¬ 
oscillation is not uncommon in electronic systems having 
negative damping due to positive feedback leading to in¬ 
stabilities. The oscillation time-period increases at a lower 
damping (see Fig. EJd) since it takes more time for magne¬ 
tization to traverse for a lower damping parameter. 

The research on single-phase multiferroic materials is 
still emerging, and the search for a room-temperature sys¬ 
tem that requires a low enough electric field for switch¬ 
ing the polarization is still underway. It will be interest¬ 
ing to incorporate the thermal fluctuations in the model 
to understand the consequence on magnetization dynam¬ 
ics mm- Also for a shape-anisotropic single-domain 
nanomagnet, the corresponding anisotropy needs to be in¬ 
cluded for detailed simulation [39]. 

Conclusions. — We have investigated the electric 
field induced magnetization switching dynamics in single¬ 
phase multiferroic materials. The dynamical system anal¬ 
ysis, contrary to steady-state analysis, revealed important 
intriguing phenomena of switching failures and oscillatory 


mode of magnetization. The key parameters that can 
shape the dynamics of magnetization are identified. The 
phenomenological magnetization damping turns out to be 
a key parameter that can prevent successful switching of 
magnetization. Hence, the present analysis puts forward 
an important step toward analyzing magnetization switch¬ 
ing dynamics between 180° symmetry equivalent states in 
the emerging multiferroic materials. Moreover, the anal¬ 
ysis identifies the oscillatory mode of magnetization that 
can act as a source of microwave signals. The present 
analysis can facilitate designing multiferroic materials for 
relevant technological purposes. 
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